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We theoretically investigate Josephson junctions with a phase shift of 7r in various proximity 
induced one-dimensional superconductor models. One of the salient experimental signatures of 
topological superconductors, namely the fractionalized 47 t periodic Josephson effect, is closely re¬ 
lated to the occurrence of a characteristic zero energy bound state in such junctions. We make a 
detailed analysis of a more general type of 7r-junctions coined “phase winding junctions” where the 
phase of the order parameter rotates by an angle 7r while its absolute value is kept finite. Such 
junctions have different properties, also from a topological viewpoint, and there are no protected 
zero energy modes. We compare the phenomenology of such junctions in topological (p-wave) and 
trivial (s-wave) superconducting wires, and briefly discuss possible experimental probes. Further¬ 
more, we propose a topological field theory that gives a minimal description of a wire with defects 
corresponding to 7r-junctions. This effective theory is a one-dimensional version of similar theories 
describing Majorana bound states in half-vortices of two-dimensional topological superconductors. 

PACS numbers: 03.65.Vf, 72.15.Nj 


I. INTRODUCTION 

In a 2001 paper, Kitaev predicted the existence of un¬ 
paired Majorana zero modes (MZM) localized at the ends 
of a proximity effect induced one dimensional (ID) p- 
wave superconductor 1 . The Bogoliubov deGennes (BdG) 
mean field Hamiltonian of this “Kitaev chain” is dis¬ 
tinguished from a trivial gapped ID system by a Z 2 - 
invariant. This topological invariant can be expressed in 
terms of the Pfaffian of the Bloch-Hamiltonian in the Ma¬ 
jorana representation. In the more recently established 
periodic table of topological states-®, this invariant is 
located in the column for dimension d = 1 in the row 
for symmetry class D , i.e. the class of superconductors 
without any additional symmetries. 

A single channel nanowire with Rashba spin orbit cou¬ 
pling, in proximity to a bulk s-wave superconductor, and 
subject to an external magnetic field, has been one pro¬ 
posal for an e xperimentally viable realization of the Ki¬ 
taev chairP 7 . A different approach taken is a magnetic 
impurity chain on top of a superconductor®®} 

Due to their charge-neutrality and non-magnetic na¬ 
ture, the unpaired MZMs are not easy to detect. The two 
main proposed signatures are a zero bias anomaly when 
the wire is coupled to a normal metal lead, and an anoma¬ 
lous 47r-periodic Josephson effect. Experimental evidence 
for the zero bias anomaly has been reported by several 
experimental groups^"®. However, it is fair to say that 
alternative explanations for robust zero bias r esonan ces, 
not related to MZMs, have also been proposecfTED. So, 
in spite of a huge experimental effort, there is still no 
uncontested experimental realization of a ID topological 
superconductor. The search for alternative observable 
signatures of this state thus remains a key challenge. 

In this paper we investigate the spectroscopy of sub¬ 


gap modes in different types of Josephson junctions in 
some detail, and ask to what extent this might pro¬ 
vide such an alternative signature. Apart from the fre¬ 
quently considered junctions, in which the order param¬ 
eter changes sign by going through zero, we also con¬ 
sider junctions for which the phase of the order parameter 
winds, while the amplitude stays constant. The sub-gap 
modes in these junctions can, at least in principle, be 
detected by standard probes sensitive to the density of 
states, and in particular scanning tunneling spectroscopy. 
Since ordinary s-wave superconductors can also have sub¬ 
gap modes in Josephson junctions, we want to identify 
spectral features that are specific to the Kitaev chain. 

We note that several other studies, complementary to 
ours, have investigated various aspects of Josephson junc¬ 
tions in topological wires 1 '"" 1 - 1 . 

The 47r-periodicity of the Josephson effect occurring 
in a junction between two Kitaev chains was pointed out 
already in Ref. Q] (see, e.g. Ref. (2Q| for a detailed dis¬ 
cussion). Closely related to this 47r-Josephson effect is a 
characteristic level crossing between two sub-gap states 
associated with a change in the fermion parity of the 
many body ground state. This level crossing is accompa¬ 
nied with a fermionic zero energy state localized in the 
junction region. Here, we study the physics of such junc¬ 
tions in both s- and p -wave paired wires from a topolog¬ 
ical perspective, focusing in particular on the nature of 
the previously mentioned (Dirac) zero mode located at a 
7r-junction. We recall how the level crossing at phase n is 
protected by an additional pseudo time reversal symme¬ 
try (PTRS) which is present in Kitaev’s minimal model 1 
for the Majorana wire if the pairing field is real (up to 
a constant phase). This additional symmetry, which is 
well known to refine the Z 2 parity to an integer winding 
number- 1 2 -( also protects the localized zero mode in the 
junction region. A major part of our present work is de- 
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voted to the study of the more general case where the 
phase of the superconducting order parameter is allowed 
to wind in the complex plane in the junction region, thus 
locally breaking the PTRS. We compare the properties of 
the 7r-junction in the topologically non-trivial p-wave case 
with those in the trivial s-wave case. Even in s-wave su¬ 
perconductors, there can still be localized sub-gap modes 
at a Josephson junction, but there is no protected zero 
mode. 

Although Kitaev’s original lattice model can be solved 
numerically for rather large systems and arbitrary junc¬ 
tion profiles, it is nevertheless interesting to verify the 
presence of the sub-gap modes, and in particular the 
zero mode, in the junction by analytical means. To 
achieve this, we linearize the spectrum around the Fermi 
points to obtain a Luttinger model, augmented with 
anomalous, charge non-conserving terms, which is es¬ 
sentially equivalent to the Su-Schrieffer-Heeger model for 
polyacetylene^l, with the Josephson junction playing the 
role of the famous domain wall soliton^. This allows 
us to find an analytical solution for the zero mode, and 
also, for a special order parameter profile, for the full 
sub-gap spectrum. Although the topological properties 
of this linearized model differ from those of the original 
Kitaev chain, we present both theoretical and numerical 
arguments for them describing the same physics. First 
we compare with an alternative linearized model (called 
below the ”V-shape model”) which does have the same 
topology as the Kitaev chain. Since this model differs 
from the first linearized model only at high (~ A) en¬ 
ergies it gives theoretical support for our claim that the 
extended Luttinger model indeed describes the low en¬ 
ergy features of the Kitaev chain. Secondly, the analyt¬ 
ical results from this model agrees extremely well with 
the numerical results obtained by directly diagonalizing 
the Kitaev chain. 

Experimentally, the most obvious way to induce a junc¬ 
tion in the wire such that the order parameter changes 
sign, is by proximity effect from a bulk superconduc¬ 
tor with a real, sign changing order parameter already 
present - this is the original scenario considered by Ki¬ 
taev. In such a junction, it is natural to assume that the 
induced order parameter in the wire remains real also 
in the junction region, and thus has to vanish at some 
point. An alternative way to introduce a junction is to 
place the wire on top of a bulk superconductor through 
which a current is flowing between two external leads 
placed below the wire. The resulting phase gradient is, 
by proximity, also present in the wire. The resulting 
’’phase winding junction” violates the PTRS, and the 
zero energy state is transformed into a finite energy sub¬ 
gap state. 

When discussing topological phases, it is interesting 
to ask what is the minimal model that will encapsulate 
the topological properties of the phase, and in particu¬ 
lar those of the elementary excitations. Important ex¬ 
amples are the Chern-Simons theories describing various 
Quantum Hall liquidfl^, and the BF theories describing 


superconductors and topological insulator^IlSSl j n thg 
present case, the elementary excitations carrying topo¬ 
logical charge are widely separated 7r-junctions at fixed 
positions, and we show that the linearized model, in the 
background of these 7r-junctions can be mapped onto 
a Dirac equation with a Goldstone-Wilczek type mass 
tempD. 

We take this as a starting point for constructing 
an effective topological field theory describing the soli- 
tons and their associated zero modes, and comment on 
similar attempts in the case of the 2D topological super¬ 
conductor. 

This article is organized as follows: In the next section 
we first define the models that we shall study. In section 
ill] we study junctions with a real order parameter for the 
different models and with both analytical and numerical 
approaches. Section |lV| contains a similar analysis for the 
phase winding junctions with constant absolute value of 
the order parameter, but in this case we have to rely 
more heavily on numerics. Section [V] briefly discusses 
possible experimental configurations to study the physics 
of topological 7r-junctions, and finally, in section |VT| we 
construct the topological field theory referred to above. 
We end with a few concluding remarks. Some technical 
points, and in particular a discussion of the rather sub¬ 
tle fc-space topology of the linearized models, are put in 
appendices. 

II. MODELS 

To set the stage for our analysis, we here first define 
the various models for the superconducting wires studied 
below. 


A. The p- wave wire 

The Hamiltonian for a spinless (or spin polarized) ID 
p -wave superconductor can be written as 

H p = J dxHp = J p,)ip 

+ A p (x)ip(-id x )ijj + A* p (x)^ (-id x )^] , (1) 

where ijj is a fermionic field (for simplicity we sometimes 
suppress the ^-dependence), p is the chemical potential 
and A p (x) = A(x)/kp is the dimensionless p -wave super¬ 
conducting order parameter. The order parameter, A(x) 
is defined such that, for constant A, the energy gap is 
2A. 

By discretizing the Hamiltonian (JT|) we get the Kitaev 
chain model 1 

JV-l 

h k = (~t( a ] a j+ 1 + a )+i a i ) 

o=o 

+ A o a o a j +1 + A j a )+i a ) - p(a]aj - j)) • ( 2 ) 
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Here, the at are (spin polarized) fermion operators, and 
we have set the lattice parameter to unity for simplicity. 
The hopping parameter is denoted by t, /i is the chemical 
potential and A j is the superconducting order parameter 
which can be position dependent. These parameters are 
related to those in the continuum model by t = 1/(2 m ), 
and [i = p — 2t. 

To write the Hk in momentum space (assuming con¬ 
stant A), we introduce the Narnbu spinor ’Pj' = (a \, a_fc), 
in terms of which, 

Hk = Y,*kUK(k)* k , 
k 

with Hk (fc) given by 

H.K(k) =(—/Li/2 — t cos(fc))r~ — Re(A) sin(fc)r y (3) 
+ Im(A) sin (k)r x: 

where the Pauli-matrices r,; act in the particle-hole spinor 
space. It is known 1 !, that for a constant order parameter, 
i.e. A j = A, the Kitaev chain resides in a topological 
phase when A/0 and |/i| < 2|t|. 

B. The s-wave wire 

As discussed in the introduction, we will compare the 
results for the topological wires with their topologically 
trivial, s-wave paired, counterparts. These trivial wires 
are described by the continuum Hamiltonian, 



+ J dx(A s (x)ip\4’l + A*(x)ipiip t )), 


where <r is the spin index, ijj^, i/j± are fermionic fields, p 
is the chemical potential and A s (x) is the position de¬ 
pendent s-wave order parameter. 

C. Two linearized models 

To capture the behavior of the above models close to 
the Fermi energy, we expand i/j into fields containing only 
low energy degrees of freedom. We consider two differ¬ 
ent ways of doing this, which give the same low-energy 
physics, but differ in their topological properties. 

1. Luttinger like model 

There is a standard way to linearize that is illustrated 
in Fig. [TJc), where the parabolic band is replaced by a 
Dirac like dispersion relation. Just as in the Luttinger 
model, we have extended the spectrum by adding un¬ 
physical ’’positron” states. In the Luttinger model, a 


gap can be opened by 2 kF processes that scatter elec¬ 
trons between the two Fermi points. In our case a gap 
is opened by charge non-conserving processes that cre¬ 
ates or destroys a Cooper pair formed by two electrons 
at different Fermi points. 

Formalizing this argument we first define, 

^=^( eikFX <P++e- ikFX <P-), 

where, Pf = y/2mp is the Fermi momentum, and <p+ 
and are right and left moving fermion fields respec¬ 
tively. Inserting this expression into ([!]) , neglecting terms 
~ e ±2lkp , we obtain 

H hin = ^ J dx(-iv F <p+d x (p + + iv F <P-d x ip- + 

+ 2(A(x)<p_<p + + A*(x)<p i + <pl)), ( 5 ) 

where the Fermi velocity is vf = Pf/tu- The quadratic 
dispersion, e(k) = ^ — p, is thus effectively replaced 
by two bands, corresponding to the right and left mov¬ 
ing linearized fermionic fields, with dispersion relations 
e±(k) = Avpk — p. In terms of the momentum q 
relative to the respective Fermi momenta, this reads 
e ±(9) = ±UFg. The superconducting order parameter 
couples these right and left moving fermions. By intro¬ 
ducing the spinor \fd = ^ + ^pp~) (the factor i is for no- 
tational convenience) we get, after integration by parts, 
the linear Hamiltonian 

H U n = J dx^Hunix)^ 

with 

H hin = ^(- iv F r z d x - 2(Re(A (x))t v - Im(A(x))r a; )), 

( 6 ) 

where the Pauli matrices now act in right-left spinor- 
space. The pairing term is taken so that the gap for 
constant A coincides with that in the previous models. 

In the following it will be important that, after rescal¬ 
ing vf by 1/2, the Hamiltonian ([b]) is identical to the 
one used by Takayama, Lin-Liu and Maki (TLM)P2, to 
describe the zero energy soliton solutions of the poly¬ 
acetylene chain model introduced by Su, Schrieffer and 
Heeger (SSH)P. We shall therefore refer to it as the TLM 
model. 

A linearized version of the trivial wire described by 0 
can be constructed in a similar fashion, but with a four 
spinor containing the left and right components of the 
two spin polarizations. For details, see Appendix [B| 

At this point we should point out that the Hamilto¬ 
nian Hun presents conceptual problems, and does not fit 
easily into the usual topological classification. The rea¬ 
sons are as follows: In Appendix A we show that as a 
consequence of the spectrum in Fig. [ljc) extending from 
plus to minus infinity, the fc-space topology is not well 
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FIG. 1: Schematic dispersion relations for the free fermion 
models corresponding putting A = 0 in (a) H p solid line, (b) 
Hv dashed line, and (c) "Hun dotted lines. 

defined. Also, the particle-hole symmetry is not a conse¬ 
quence of a redundancy due to an artificial doubling of 
a band. Rather it follows from extending the linear dis¬ 
persion to arbitrary large negative energies. If we were 
to add band bending corrections to this model we would 
break the particle-symmetry which again would change 
the topological classification of the model. This situa¬ 
tion is unsatisfactory since it raises questions about the 
validity of linear approximations, and in particular the 
use of the TLM model, for analyzing the Kitaev chain. 
To resolve this we shall now present an alternative model 
that resolves the problems related to topology and dou¬ 
bling, while retaining a linear spectrum. Having shown 
the existence of such a model, we can safely continue to 
use i?Lin br the subsequent discussion. 

2. V-shape model 

First we replace the parabolic band in Fig. |TJa) with 
a V-shaped band, with dispersion e v (k) = \k\vF — jd, as 
shown in Fig. |TJb) . 

Next we express the full field i/i(k ) in terms of the low 
momentum fermion fields tp±{k) 

m = ^(e ik r*v + (k)m+e- lkFX V-(m-k)), ( 7 ) 

where 'd(fc) is the step function. In order to write a BdG 
Hamiltonian, we first define 

X(k) = -^(e ik r*<p + (k)d(-k) + e- ik '*<p-(k)m ), (8) 

and the Nambu spinor Next we substi¬ 

tute |7| and ([8]) in the expression for H p , and disregard 
the rapidly oscillating terms ~ to get 

Hv(fc) = [(-fi + v F \k\)r z + A sgn(fc)r y ] <f> , (9) 

where again the Pauli matrices t, act in the Nambu space. 
As usual, this amounts to a doubling of the spectrum, 
and this redundancy is manifested in the particle-hole 
symmetry of "H v which cannot be broken. The pairing 


term ~ A (which is assumed to be real) is such that it 
gives rise to the same gap as the original Hamiltonian H p 
for constant A. 

By inspection, we see that the dispersion relation for 
7~L v (k) has an unphysical 2A jump at k = 0. This dis¬ 
continuity can be regularized by smoothening the tip of 
the V-shaped band, and this will in fact be necessary 
when we analyze the topological properties in Appendix 
A. Such a regularization will however necessarily yield a 
more complicated model, that is only amenable to nu¬ 
merical solutions, in spite of having a very simple low 
energy limit. We will not pursue this since, this model is 
of interest only to demonstrate the existence of a consis¬ 
tent model with a linear spectrum, and good topological 
properties. 


III. JUNCTIONS AND SOLITONS 
A. Topological properties 

We start our discussion of 7 r-junctions in ID super¬ 
conductors, by reminding the reader about which differ¬ 
ent topological superconductors are possible in ID sys¬ 
tems. To do this, we recall the topological classification 
of non-interacting fermion systems 2 4 , where the possible 
topological phases are classified according to their non¬ 
unitary symmetries, viz. time-reversal symmetry (TRS) 
T and particle-hole symmetry (PHS) C (we note that 
the PHS is technically a spectral constraint rather than 
a physical symmetry. However, we here chose to follow 
the widely adopted terminology of Ref. n. 

In this paper, we consider superconductors in one di¬ 
mension without spin rotation symmetry. The BdG 
structure of the Hamiltonian entails a built in algebraic 
constraint rooted in the fermionic algebra of the field 
operators that can formally be viewed as a PHS with 
C 2 = +1. In the absence of time-reversal symmetry, 
i.e. for class D, the superconductor is either topolog¬ 
ically trivial, or non-trivial, depending on the value of 
the Z 2 invariant. In the latter case the wire supports 
MZMs at both endsP. In the case of time-reversal sym¬ 
metric superconductors, with T 2 = —1, i.e. in class Dill, 
the situation is similar, but in this case, the topological 
phase exhibits a Kramers-degenerate pair of MZMs at 
both edges, see, e.g., Refs. l30|31l ) Finally, if the sys¬ 
tem respects the PTRS T 2 = +1, i.e. for class BDI, the 
different topological phases are distinguished by an in¬ 
teger winding number, giving an infinite set of different 
topological non-trivial phases. 

The p -wave wire, ([lj or ([ 2 ]), will in general, i.e. when 
we allow both the hopping and the order parameter to 
be complex, belong to symmetry class D, which means 
that it can either be in a trivial phase, or in a topological 
phase. In the lattice model, the former happens for \p\ > 
2\t\, while the latter occurs for \p\ < 2\t\, with |A| 7 ^ 0. 

If both t and A are real, the Hamiltonian © is also 
pseudo time-reversal symmetric (here T is simply com- 
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plex conjugation, so trivially T 2 = 1), and in this case, 
the possible topological phases are labeled by an integer, 
corresponding to a winding number (see Appendix 
Kitaev’s model with a constant order parameter ex¬ 
hibits three of these phases, namely the trivial one (when 
|/r| > 2|f|), as well as two non-trivial ones, both occurring 
for \fi\ < 2\t\, one with A > 0, the other with A < 0. 

This means that for real t and A, the Kitaev chain 
can harbor an interesting junction, by allowing the order 
parameter to change from — A to +A in a finite region, 
corresponding to a 7r-junction. Just as the edge of a 
Majorana wire hosts a MZM, because it constitutes the 
boundary between a topological phase and the trivial vac¬ 
uum, the 7r-junction we consider here will also support 
zero modes. Since the difference in winding number be¬ 
tween the two neighboring topological phases is two, we 
expect twice as many zero modes in comparison to the 
edge of the Kitaev chain. Below we show that this is 
indeed the case, irrespective of the precise ^-dependence 
of the order parameter. 

We already mentioned the problems related to properly 
define the fc-space topology for the TLM model, and how 
they are resolved by an alternative linearization scheme. 
The details are given in Appendix A, but we should here 
again stress that the outcome of this analysis is that we 
can safely use the TLM model to discuss the topological 
properties of the Kitaev chain. 

B. The 7r-junction as a soliton 

Although the k space argument for topology of the lin¬ 
earized model Hun given in Appendix[A]is compelling, it 
is important to find out how well the TLM model ([b]) re¬ 
ally captures the topological properties of the full model 
([2]). To do this, we first briefly recall how a Dirac Zero 
Mode (DZM) arises in the TLM modeP^l, and then com¬ 
pare it with the zero mode arising in the full model ([ 2 ]), 
in the presence of a junction, at which the real order 
parameter A changes sign. 

The presence of the DZM in the case of a real order pa¬ 
rameter is most easily demonstrated in the TLM model, 
and from TLun we get the BdG equations 

-(— ivpd x u{x ) + 2iA*(x)v*(x)) = cu(x) 

-(ivpd x v*{x) — 2iA(x)u(x)) = ev*(x) . (10) 

For real A (a;), and taking e = 0 since we are interested in 
the zero modes, these equations are easily decoupled by 
introducing f±(x) = u(x) ± v*{x). For a 7r-junction that 
interpolates between a negative constant A_ for x <C 0 to 
a positive constant A + for x 0, one finds the solution 
f+(x) = NeA(x )dx — Q Here, we will 

consider the special profile A(x) = A 0 tanh(a;/£), that 
gives rise to the analytical solution— 

/+ (z) = N 0 sech(x /£) c/?0 /_ (x) = 0 , 


t 

Ao 



MLS Error 

10.0 

1.0 

0.0 

0.999418 

2.2002 • 

10~ 8 

8.0 

1.0 

0.0 

0.999098 

6.7994 • 

10 -8 

5.0 

1.0 

0.0 

0.997777 

7.5333 • 

10- 7 

2.0 

1.0 

0.0 

0.991853 

1.2267- 

10” 4 

1.0 

1.0 

0.0 

1.071428 

6.8020 ■ 

10~ 3 


TABLE I: The fitting parameter A, compared to its analytic 
value in the TLM model £o> as well as the error of the fit (using 
the Method of Least Squares), for a system with N = 1001 
sites and various values of t. 


where £0 = vf/( 2Ao) and No a normalization constant. 

We compare the DZM of the TLM model to the 
full model, by considering the discretized version Aj = 
A 0 tanh(j/£) of the TLM profile A(x) = A 0 tanh(x/£) 
in ([ 2 ]). By choosing the width £ not too large, this de¬ 
termines the order parameter to — Aq at one end of the 
chain and +Ao at the other end, generating a domain 
wall (between sectors with different winding numbers) at 
the center of the chain. We set the junction parameter 
£ = £ 0 , and fit the resulting zero mode to the TLM so¬ 
lution g(x) = sech(x/A), with A used as a fitting param¬ 
eter. In Tab. [I] we show some representative results and 
in Fig. [2] we display a typical result for the probability 
density of the DZM located in the junction, as obtained 
from the Kitaev chain. Evidently, the TLM model cap¬ 
tures the the properties of the DZM in the junction region 
of the Kitaev chain very well. 



Site number 


FIG. 2: The probability distribution (black dots) of the zero 
energy mode located in the junction region of the Kitaev chain 
with order parameter profile ~ tanh(x/£o). The parameters 
used are: t = 10.0, Ao = 1.0 and /.t = 0.0 yielding £0 = 10. 
The fit (red, dashed line) is made by the method of least 
squares, and resulted in A = 9.994, in good agreement with 
the value of £o- The number of sites is N = 1001; the figure 
only shows the central region. 


( 11 ) 
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IV. PHASE WINDING JUNCTIONS 


The piecewise solutions are given by 


In this section, we extend the previous discussion to 
the case of junctions with a complex order parameter 
where the phase winds in a finite segment of the wire. 
For simplicity we shall assume that the absolute value 
|A| remains constant. In this case we can find an analyt¬ 
ical solution in the linearized model by taking a simple 
winding profile, while our numerical analysis easily gen¬ 
eralizes to more general profiles. 

Although a complex order parameter breaks the PTRS 
and with that the chiral symmetry that protects the 
DZMs in the junctions with real profiles, one would still 
expect that the low energy theory should not differenti¬ 
ate a rapid winding of the phase from 0 to 7r from a sharp 
step in the magnitude of A. We now demonstrate that 
this intuition is correct, and that low energy modes per¬ 
sist even in the case of phase winding junctions. Again 
it is easiest to start from the linearized model. 


! aie ra x < —a 

(a^e KX + a 2 e~ KX ) jo) < a (17) 
tt 3 e _KI x > a 


where n = y C 0 2 — e 2 and k = y C 0 2 — (e + 7r/(4a)) 2 . 
To obtain a normalizable solution, we must take n < 
0, or |e| < 1/Co, implying that the (sub-gap) solution 
is localized in the junction region. From the matching 
conditions for the wave function and its derivative, one 
can infer that there is no solution when k is real. An 
imaginary k requires that e > 1/Co — 7r/(4a), so localized 
sub-gap modes are possible in the energy range 1/Co — 
7r/(4a) < e < 1/Co if a > £ 0 ^/ 8 , or in the whole gap 
region -1/Co < e < 1/Co if a < CoV 8 - 

For imaginary k, the matching conditions have a solu¬ 
tion if the following constraint is satisfied 


A. sub-gap states in the linearized model 

We consider an order parameter with an ^-dependent 
phase 


A = A o e I0(x) 0(-oo) = 0 0(oo) = fir, (12) 



where Ao is a positive constant, 9(x) is continuous and 
/ is some real number. The BdG equations (10) then 
become, 


— id x u(x) + i/Coe ie ^v*(x) = eu(x) (13) 
id x v*(x ) — i/t;oe l6 ^u(x) = ev*(x ), 

where Co = v f/( 2A 0 ) and e = 2 e/vp. From the first 
equation we have 

v* ( x) = -i£ i0 e l9{ ' x) (e + id x )u(x) (14) 

and substituting this into the second, we get 

[d x + i{d x 0)d x + e(d x 6 ) + (e 2 - Cc/ 2 )] u(x) = 0 . (15) 

This equation cannot be solved analytically for a general 
profile, but for the case of 


9 k (x) = 


0 x < —a 

1*1 < a 

fir x > a 


(16) 


we can solve (15) in the three regions and then match 


the solutions. Just as in an ID Schrodinger problem in 
a piece-wise constant potential, this is done by matching 
the function and its (logarithmic) derivative. We focus 
on the case / = 1, which corresponds to a 7r-junction 
where A changes sign, but the analysis below can easily 
be extended to junctions with arbitrary phase winding. 


Upon analyzing this equation, one finds that even for ar¬ 
bitrary small a, there is always at least one solution. The 
energy of the associated bound state is always positive, 
but approaches zero in the limit of small a. Upon in¬ 
creasing a, more and more bound state solutions appear. 
In order to have at least p + 1 bound states, a should 
satisfy a > 

Before turning to the numerical results, we briefly dis¬ 
cuss the case of general phase winding, ie., we allow / 
in © to be arbitrary. For / arbitrary small, one finds a 
bound state, with an energy slightly below the band gap, 
e < 1/Co- Upon increasing /, the energy of this bound 
state decreases towards e = — 1/Co- In the mean time, 
more bound states appear at the gap edge e = 1/Co- In 
the limit of large /, the energies of the bound states be¬ 
come periodic in /, with a period of 2, i.e., a period of 
2-7T in the winding angle. Finally, in the limit of a very 
short junction, we find that for / an odd integer, there 
is a bound state at e ss 0, while for / an even integer, 
there are two bound states with energy e « ±l/Co- In 
the former case, the junction behaves as a 1 junction 
with a real order parameter, while the second case is 
equivalent to not having a junction at all. This is con¬ 
sistent with the topological discussion above, although 
we should point out that there are no topological reasons 
why the phase junction should behave as a real junction 
in the short junction limit. We next compare some of the 
results of this section with numerical simulations in the 
Kitaev chain and in the full s-wave model. 
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Analytic sol. 

Linear model 

Full model 

0.940199 

0.940201 

0.940156 

0.956549 

0.956556 

0.956385 

0.981316 

0.981324 

0.981002 


TABLE II: The energies of the first three bound states in a 
p-wave 7r phase winding junction. The parameters used for 
these calculations are: t = 10.0, A = 1.0, p = 0.0, a = 120, 
N = 800. Analytical, linear model and full model refer to 
the equations ( fl8] ), ([To]) and respectively. Note that the 
linear model values are just a measure of how well analytic 
solution describes the discretized linear model, while the full 
model values describe how well the linearization captures the 
low energy degrees of freedom. 



B. Comparison with the Kitaev chain 


Starting with the Kitaev chain model given in (pi), we 

(j+a) 1 —' 

take the profile A j = \A 0 \e X7r ^~, so that over a seg¬ 
ment of length 2 a, the phase increases linearly from 0 to 
7r. Effectively, this amounts to changing the sign of A just 
as in the previous section. Using this profile, we numeri¬ 
cally calculated the energy of the low lying fermion states 
both for the Kitaev chain and the linearized model, using 
a range of parameters. Typical results are shown in Tab. 
PI where the agreement between the first two columns 
is a measure of the precision of our numerics, and the 
good agreement with the third column again confirms 
that the linearized model faithfully describes the full Ki¬ 
taev chain. We have also compared the numerical wave 
functions for the low lying states in the Kitaev chain, 
with the analytical expressions ( |17[ ) and again found ex¬ 
cellent agreement. 

Next we studied what happens when the length of 
the phase winding 7r-junction shrinks. In Fig. |3j which 
shows our result for the p -wave case, we see clearly how 
a state that is close to the gap for large junctions comes 
down, and becomes a zero mode for the shortest junc¬ 
tions (which essentially amounts to a sign change between 
two lattice points). This supports the heuristic argu¬ 
ment, given earlier, that a short phase winding 7r-junction 
should have properties very similar to the one where A 
remains real but changes sign. The corresponding s-wave 
setup is depicted in Fig|4j where no zero modes need to 
be formed in the short junction limit. 


These results give additional confirmation that the low 
energy properties of junctions made by Kitaev chains can 
be captured by the linearized model in ([5]), and in the 
Section VI we construct a topological field theory, which 
captures the same physics. 


FIG. 3: The energy spectrum of the Kitaev chain as a function 
of its 7r phase junction length (in units of the lattice param¬ 
eter).. A junction length of 1 means that the phase jumps 
from 0 to 7r from one site to another. Note that the zero 
energy states that represent MZMs located at the end points 
of the chain have been omitted. In addition, two new zero 
modes are formed as the junction length shrinks, effectively 
imitating a real 7r-junction. The parameters used are t = 2.0, 
Ao = 1.0, p = 0.0 and N = 200. The spectrum is displayed 
in a low energy regime. 



FIG. 4: The energy spectrum of the s-wave wire as a function 
of the length of its 7r phase junction (in units of the lattice 
parameter). A junction length of 1 means that the phase 
jumps from 0 to n from one site to another. Because the s- 
wave wire is topologically trivial, no zero modes form, even in 
the limit of a short junction. The parameters used are t = 2.0, 
Ao = 1.0, p = 0.0 and N = 200. The spectrum is displayed 
in a low energy regime. 


V. HOW TO EXPERIMENTALLY PROBE 
TOPOLOGY BY A tt-JUNCTION 

Most of the experimental effort in studying the topo¬ 
logical wires has been aimed at detecting the MZM at the 
edges. But as mentioned in the introduction, the pro¬ 
posed signatures for these modes can also be emulated 
by other effects. It is thus interesting to consider other 
signatures for the wire being in the topological phase, 






































































































and here we suggest the possibility of using the DZMs at 
7 r-junctions as such a probe. For this idea to be useful, 
we not only need a way to experimentally realize such a 
junction and detect the associated fermionic zero modes, 
but also a clear signature for the topological phase. We 
shall consider both junctions with topologically protected 
zero modes, and phase winding junctions. We begin with 
the latter. 


A. Phase winding 7r-junction 

One way to make a phase winding junction is to put 
a wire of the type used in previous experiments on top 
of a s-wave superconductor through which a current is 
driven between external leads placed close to the wire. 
By the the relation V(f> ~ J, where <j> is the supercon¬ 
ducting phase and J is the current, one can arrange for 
a 7 r phase difference between the leads, which will, by 
proximity, be imprinted on the wire. An experimentally 
more challenging task is to probe the fermion spectrum 
at the junction. An obvious possibility is to use a tun¬ 
neling contact weakly coupled to the wire, or a scanning 
tunneling microscope. 

From the previous section it would appear that a good 
signature for the p-wave pairing phase would be the pres¬ 
ence of an almost zero mode in the junction region. Un¬ 
fortunately, the situation is not very clear since an s-wave 
pairing would have a similar signature. Fig. [4] is similar 
to Fig. [3j but for s-wave pairing. Also here we find a 
low-lying sub-gap state for short junctions, and although 
it does not come all the way to zero, it is not clear that it 
could be distinguished from the p -wave case. Clearly one 
would need much more detailed studies of more realistic 
microscopic models in order to resolve this question. 


B. Real 7r-j unction 

As already pointed out, in a 7r-junction with a real 
order parameter (that must go through zero) the zero 
energy Dirac mode is always present when the supercon¬ 
ductor is in the topological phase. For the trivial s-wave 
case, there is no such protected zero mode, but the spec¬ 
trum of the subgap modes does depend on the profile of 
the order parameter at the junction (and on the other pa¬ 
rameters, such as the chemical potential). Importantly, 
there can be junction modes with zero energy, that can 
be described by the TLM model we studied above, for 
certain choices of parameters. For example, putting the 
chemical potential in the band middle (p = 0 or equiva¬ 
lently p = 2t as measured from the bottom of the band), 
there are localized modes with zero energy, regardless of 
the junction length. But these states can be gapped out 
in the short junction limit by lowering the chemical po¬ 
tential to the vicinity of the band bottom. This feature 
is demonstrated in Fig. [5] and contrasted with the cor¬ 
responding p-wave system in Fig. [g] In the latter case, 


the topology protects the zero mode, regardless of the 
junction length, as long as the chemical potential lies in 
the band, so that the system is in the topological phase. 

We note that for wide p-wave junctions, £ > 30, there 
are additional subgap modes with finite energy which are 
not in the range of ^-values in Fig. [6] 



Junction width 


FIG. 5: The energy spectrum (bulk states in black and subgap 
states in red) of the s-wave wire (a discrete version of Eq. 0 ) 
as a function of the width of its real 7r-junction (in units of 
the lattice parameter). Modes with zero energy exist only 
in the limit of a wide junction, and are gapped out in the 
short junction region due to the low chemical potential. The 
parameters used are t = 1.0, Ao = 1.0, p = 1.9 (p = 0.1) and 
N = 200. The spectrum is displayed in a low energy regime. 



FIG. 6: The energy spectrum (bulk states in black and subgap 
states in red) of the p-wave wire (Eq. 0 ) as a function of the 
width of its real 7r-junction (in units of the lattice parameter). 
The Dirac zero mode is topologically protected and exists for 
short and long junctions since the wire is in its topological 
phase. Subgap modes with finite energy are not present in 
the junction length regime displayed here. The parameters 
used are t = 1.0, Ao = 1.0, p = 1.9 (p = 0.1) and N = 200. 
The spectrum is displayed in a low energy regime. 
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There are at least two possible ways to experimentally 
realize a junction of this type. The most direct would be 
to implement a 7r-junction in the underlying s-wave su¬ 
perconductor, but the problem here is that it is not easy 
to determine the A profile in the junction. An alterna¬ 
tive way is to put the wire as a bridge connecting two 
different s-wave superconductors that are held at differ¬ 
ent values of 6 , for instance by a SQUID geometry. In 
this case one might calculate the A profile by a realistic 
modeling of the wire. 

Clearly both options need to be studied in more de¬ 
tail before any definite conclusion can be made about 
the feasibility of using sub-gap junction spectroscopy for 
probing the topological nature of the wire. 


VI. A TOPOLOGICAL FIELD THEORY FOR tt 
PHASE WINDING JUNCTIONS 


As in the previous sections we consider A as given by 
the background s-wave superconductor by proximity. We 
generalize the previous discussion somewhat by consider¬ 
ing a wire with many widely separated 7r-junctions of the 
phase winding type. Since the bound states are exponen¬ 
tially localized, such configurations will support subgap 
modes at each junction. We stress that we consider a 
fixed junction configuration, given by the background s- 
wave order parameter, and derive an effective theory for 
the fermions. We can, using the same formalism, also de¬ 
scribe adiabatic motion of the junctions, but they cannot 
be considered as bona fide itinerant particles. 

Our starting point is the Lagrangian formulation of 
the linearized version of the p-wave superconductor ©. 
which is given by 

C = if(i$- g(x)e M{x) ^ if , (19) 


where the functions d and g are related to the supercon¬ 
ducting order parameter by A = ge' ( !? + 7 ' 7 2 ) ) and where 
the Dirac matrices are related to the Pauli matrices by 
( 7 o> 7 x) = (cU, —ia y ) so 75 = er 2 . We have in our deriva¬ 
tion relabeled the spinor to if, in accordance with stan¬ 
dard notation. Furthermore, we have put vf = 2, consis¬ 
tent with the linearization of ([l]). Note that for real A, 
the energy gap A is nothing but the mass in the Dirac 
equation. 

The aim here is first to derive a bosonic form of the 
Lagrangian (19), and then to extract an effective action 
that describes the physics of the bound states on the soli- 
tons. This theory is topological in the sense that it does 
not have any bulk degrees of freedom, but only describes 
the quantum mechanics of the bound states residing on 
the solitons. Before embarking on this exercise, we will 
put it in context, and view it as part of the more chal¬ 
lenging problem to formulate topological theories in the 
presence of fermionic zero modes. 

As mentioned in the introduction, the Kitaev chain is 
a one-dimensional cousin of the two-dimensional (2D) p- 
wave superconductor, and the Majorana states located 


on the interfaces between the normal and topological 
phase of the model can, by employing geometries with 
junctions, be used for quantum computing. In this con¬ 
text the quantum mechanics of the zero modes is clearly 
very interesting. In the 2D case, there are two candidates 
for a topological field theory that describes the braid¬ 
ing of vortices with Majorana zero modes. One is based 
on a SU(2) Chern-Simons theorjEH, while the other em¬ 
ploys an abelian BF theory coupled to a single Majorana 
field 33 . In this 2D case, the vortices are in principle itin¬ 
erant, but are in practice often pinned to impurities. In 
this latter case there is a close analogy with our system 
of fixed, or adiabatically moving, 7r-junctions. 

None of these effective theories just mentioned has 
been derived from a microscopic description, but are ob¬ 
tained from general principles based on symmetry and 
scaling. The effective topological theory for the fermionic 
bound states on solitons that we shall describe shortly, is 
closely related to the second of the 2D topological theo¬ 
ries that we just mentioned. An obvious, and important, 
difference is that the fermionic modes on the solitons 
are of Dirac type, so, even if fine tuned to zero energy, 
they can not be used for topological quantum comput¬ 
ing. The advantage with the present case is that it is 
more amenable to analytical treatment. Still we have 
not managed to obtain the topological theory directly 
from the microscopic model in a controlled fashion. The 
derivation presented below is therefore phenomenological 
and again based on symmetry considerations and scaling 
arguments. In Appendix 0 we do offer a microscopic 
derivation which however involves several unproven, and 
admittedly questionable, assumptions. 


A. Symmetries 


We now discuss the symmetries of (19). From this 
Lagrangian we can immediately get the vector and axial 
charge densities, 


p v = iftif = <p' + <p+ - ( 20 ) 

Pa = if^5if = ¥+¥>+ + <P-V- (21) 


which shows that the electric charge Q em = Qa = 
/ dx pa(x) in the superconductor is given by the axial 
charge (21) in the Dirac theory (19), and is thus not 
conserved, as appropriate for a superconductor. Note, 
however that (19) is invariant under the combined global 
transformation 


if -A e i/ 37 E 7 f ( 22 ) 

e^d-2/3. 


In a BdG description this corresponds to a simultaneous 
global phase change of the electron field and the super¬ 
conducting condensate (ifif). Also note that the trans¬ 
formation, 

if -a e mi5 if = —if (23) 
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is indeed a symmetry. As expected, this is a manifesta¬ 
tion of the the conservation of electric charge modulo two, 
whic h is most easily seen by noting that the transforma¬ 
tion (231 leaves the pairing terms tp^ + tp\_ and tp + tp_ in¬ 
variant. It will be important later that the vector charge 
Qv = f dx pv{x) in the Dirac theory is indeed conserved. 
Physically this is a consequence of the Cooper pairs hav¬ 
ing zero momentum, so adding or subtracting a pair will 
give identical changes at the two Fermi points^. In the 
following we shall give a bosonized version of the theory 
where it will be important to keep the correct symmetry 
pattern. 


B. Bosonization 


It will be advantageous to rewrite (19) in bosonic vari¬ 


ables using the bosonization “translation table”, (see for 
instance Ref. USD, 


ipijj ■ 
$ *75^ ■ 


K COS ip 

k simp 


(24) 


where the dimension-full parameter k depends on the 
short distance cutoff, and the scalar field <p is normal¬ 
ized so that the bosonic version of (19) is 


c = - gcositp-'d) 


(25) 


where we have rescaled g with k. The minima of the 
potential are at 


tp n = $ + 7T + n2ir = 9 — — + n2n , 


(26) 


so for large g, 9 = 9 — 7 t/2 will make small fluctuations 
around one of these (equivalent) minima. In particular, 
if 9 winds, then tp follows. From the work of Jackiw and 
RebbP^, and Goldstone and WilczelP^, we know that 
windings in the scalar field tp will describe solitons carry¬ 
ing (in general fractional) fermion number. For simplic¬ 
ity we neglect 27 t windings, and taking n = 0 in (26) we 
define the kink current as 


% = ’ 


(27) 


so the charge of the soliton that interpolates 0(x) from 
tp L to tp R is given by, 


q s = ~ ■ 


(28) 


It follows that the 7 r-junctions we discussed earlier carry 
a half unit of fermion number. 

Next we shift the field tp by tp = tp 0 + (j> to get 


C ’ = b' rdvdb »- 

+ ^(<V£) 2 cos(<^>) + ) 2 j£e llI/ d v $+ 


1 

8n 


(d,0f 


where bn is a multiplier held that imposes the condition 


(27). Since the 0-field is massive, it can be integrated, to 


yield the truly trivial topological theory, 
£eb = ^-e^dvBbu, - b^. 

Z7T 


(30) 


C. Retaining the fermion bound states 

The topological theory we just derived is however not 
always a good description of the low energy physics. This 
is most easily seen by considering the special case where 
the topological current describes widely separated nar¬ 
row 7 r-solitons. We learned in section IIV Al that these 
can support low energy fermionic bound states with en- 
ergy eo < e < A inside the gap. Since we furthermore can 
fine tune so one of these modes occurs arbitrarily close to 


zero energy, the topological theory (30) can clearly not be 
universally correct. Moving away from the eo point, but 
still having the bound state far below the bulk gap, i.e. 
eo <C g , it would still be desirable to have a theory that 
describe these low-lying excitations. What went wrong 
in the derivation of (30) is that while the bosonic fluc¬ 


tuations with energy > g were integrated, the more im¬ 
portant fluctuations changing the fermion number were 
not taken into account. We will now remedy this and 
present a model that properly includes the dynamics of 
the low-lying fermionic bound states. 

We shall first construct a model in the limit of widely 
separated point-like kinks. Any real function A that in¬ 
terpolates between ±|A 0 | at x = ±oo supports a zero 
mode. The kink, |A|?/(a;), where g is the step function, 
can be thought of as a limit of such functions, and thus 
supports a zero mode. Also, as discussed above, we get 
an approximate zero mode for constant |A| = m, and a 
rapid winding of the phase 9 an odd number of 7 r. In 
both these cases the topological current related to the 
kink can be be described by 


N 


jk(x, t) = ^ 6{X - X a ) (1, X a ) 


(31) 


where we allowed for the kink at position x a to move 
with velocity x a . 

It is now straightforward to write a Lagrangian for the 
bound states residing on the kinks, 


N 




a =1 
N 


(32) 


= ^2C(^ X a(t))i(d t - X a (t)d x )€a(t,X a (t)) 


a =1 


= J dxj^^(x,t)id^(x,t). 


(29) Combining this with the term ( |30| ), yields 

C ieb = ^d v 9 [&„ + tUdvt] - eo?Z - bpft (33) 
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where we also introduced a chemical potential eo that 
fixes the energy of the bound state. We shall take £ 
to be a complex fermionic field (otherwise it would not 
describe a single bound state), but note that it differs 
from a conventional Dirac fermion in being dimensionless. 

The first term in Lagrangian (33) is closely related to 
the topological Lagrangian for a spin-less 2D chiral su¬ 
perconductor given in Ref. l33l The main difference is 
that in the 2D case the Dirac fermion £(x, t) is replaced 
by a Majorana field 7 (x,t). In the present setting, that 
would be appropriate for a domain wall between a triv¬ 
ial and non-trivial phase of the wire. The second term 
~ eo is not topological and is present only for a com¬ 
plex field. Note that the kinetic term ~ in (33) 

has support only where the topological charge does not 
vanish, and thus there are no bulk degrees of freedom. 
The above analysis is, however, valid only for point like 
sources. The generalization to extended sources, that is 
the finite size kinks considered in the previous sections, 
is our next task. 


D. Fermion bound states in extended kinks 


Since for a static kink, the Hamiltonian in (33) is only a 


chemical potential, it can not describe the fermion modes 
on an extended kink, but instead gives a continuum of 
states at energy eo- To get a realistic low energy the¬ 
ory we must thus introduce more terms in the effective 
Hamiltonian. Following the usual logic of effective theo¬ 
ries we shall retain the lowest derivative terms that en¬ 
sure the correct symmetries. The crucial symmetry here 
is the broken global 17(1) symmetry related to the electric 
charge. In the linearized theory (19) this is the (global) 
chiral symmetry (22). Clearly terms like ^d 2 ( etc., 
are allowed, but also pairing terms like e l6 t^d x £ etc.. 
In fact it is necessary to include a pairing term in order 
to get the appropriate symmetry breaking. Putting the 
chemical potential eo to zero, the simplest possible action 
for an extended kink is, 


C& b = [ b v + b^ (34) 


with 


?*£ = (M 2 - d 2 x )t + ^ + e~ ie !?id x (?] , 

(35) 

where the mass parameter M and the pairing strength 5, 
are phenomenological parameters. 

We can simplify this Hamiltonian by performing a ro¬ 
tation of the fields: 

£ -»• e~ l9/2 Z ft -»■ e l9/2 t} ■ (36) 

This will transform the Hamiltonian ([35]) to 77^ = 


with 


«£ = 


AT 


- (| d x - | d x ef 
8Mid x 


—AT 2 + 


5Mid x 

(d x + § d x ey- 


(37) 


Next, we expand the quantum field as £(x,t) = 


En(e- iEt ° 


i{x)c 


e iEt vz 


(x)c n ), which yields the fol¬ 


lowing BdG equations for the eigenfunctions u(x) and 
v{x), 

d x + ^d x 9) 2 + Ed x 9 — M 2 ^ju(x) + 5Mid x v* (x) = 0 

(38) 

((d x - l -d x 9) 2 - Ed x 9 - M 2 y*{x) - SMid x u(x) = 0 

In the limit <5 —» 0 and under the assumption that 9 
varies slowly (i.e, we assume d 2 9 and (d x 9) 2 to be small) 
we obtain the following equation for u(x) 


[<% + (d x 9)id x + E(d x 9) - M 2 } u{x) = 0, 


(39) 


which is © in the limit where the energy E is small 
compared to M. As expected there is no continuous com¬ 
ponent in the spectrum, and the low energy part of the 
spectrum compares well with the full model with suitable 
adjustment of the model parameters. In particular, we 
should set M 2 = = 4A 2 /v 2 F , E = e = 2 c/vf- The 

requirement that E <C M then translates to e <C A, that 
is, for energies well below the gap, which is consistent 
with a zero energy bound state. 


To actually derive the effective Lagrangian (33) one 


should integrate out the high energy modes. This would 
not only give expressions for the effective parameters, but 
also provide an ultraviolet cutoff that would define the re¬ 
gion of validity of the effective model. We have not been 
able to do this in a controlled way, but in Appendix [C| it 
is shown, by manipulating path integral expressions, how 
the crucial kinetic term \e^ v d v 9^id^ can arise from the 
microscopic description. 

Finally we note that the extension of the topological 
theory (33) to the model Lagrangian (34) for the sub¬ 


gap regime, is reminiscent of the extension, proposed in 
Ref. [37; of the 2D topological theory in Ref.[33l In both 
cases the models are constructed using phenomenologi¬ 
cal and heuristic arguments, and it remains a theoretical 
challenge to find general methods to describe localized 
fermionic zero modes in the general context of topologi¬ 
cal field theory. 


VII. CONCLUDING REMARKS 

In this paper we studied several models for trivial 
and topological superconducting wires in one dimension. 
More specifically, we investigated the properties of 7 r- 
junctions, and in particular those where the phase of the 










12 


order parameter winds an angle 7 r over the junction, cor¬ 
responding to a system in symmetry class D. For this 
more general case, we find that there is no topologically 
protected zero energy mode associated with a 7r-junction. 
Rather, local breaking of the PTRS by means of the com¬ 
plex winding of the order parameter can shift the energy 
of the bound state in the junction region away from zero 
energy. This symmetry breaking is not allowed in class 
BDI, where, as a consequence, the bound state is topo¬ 
logically pinned to zero energy. We demonstrated that 
the low energy bound states in some specific cases can be 
obtained analytically and showed that these results agree 
well with numerical calculations. Most importantly, we 
discussed how our results might be used to obtain a bulk 
probe - in contrast to the common method of probing 
the edges - to distinguish a topological wire from a triv¬ 
ial one, and suggested some experimental approaches to 
this end. Finally we constructed a low energy field theory 
with a topological term describing itinerant 7r-junctions, 
and discussed its relation to theories in two dimensions. 
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Appendix A: Topology of the p-wave superconductor 
modes. 

Here, we discuss the topological properties of the vari¬ 
ous models we consider in this paper. To set the scene, we 
start by recalling the topological properties of the Kitaev 
chain, see Ref. [1]. 

Consider the model © and assume that t and A are 
both real, so that the Hamiltonian belongs to symmetry 
class BDI. The topological invariant takes the form of a 
winding number 2 !, and to show this in the present case, 
we write the fc-space Hamiltonian © as 

HK{k) = d{k ) • t , (Al) 

with t = (t x , T y . t z ). For models in class BDI, one can 
choose a basis such that one of the components of the 
vector d is zero, say d x = 0. The energy is given by 
e(k) = ±|d(fc)|, which means that for a gapped sys¬ 
tem, we have > 0. Hence, the winding number v 

around the origin of the curve in (r y , r z )-space (i.e., the 
space of Hamiltonians) swept out by d(k) as k sweeps 
through the full Brillouin zone is well defined. This 
winding number is the topological invariant character¬ 
izing the different phases. For the Kitaev chain we have 
d(k) = (0, — Asin(fc), — g/2 — t cos(fc)), and in Fig. [ 7 } we 
(schematically) show the curve d(k) in the trivial phase, 
with winding v = 0, and the two different topological 
phases, with winding v = ±1. 



FIG. 7: Winding numbers v of d(k) for the full Kitaev chain, 
in (a) trivial phase with v = 0, for 0 < t < g/2, A > 0, (b) 
topological phase with v = 1 for g = 0, 0 < t = A and (c) 
topological phase with v = — 1 for g = 0, 0 < t = —A. The 
arrows denote the direction in which k increases. 


Next we turn to the linear model Jiun hr Eq. ([ 5 ]). As¬ 
suming that A is real and constant, the momentum space 
version of the Hamiltonian ([6]) is again of the form 
with d(k) = (0, —2A, t»irfc). 

Since the fc-space is not compact, it is possible that the 
curve swept out by d{k ) = d(k)/\d{k)\ (the normalization 
is needed to obtain finite limits and is valid as long as the 
Hamiltonian is gapped) as k goes from —00 to 00 is not 
closed. This is indeed what we find in Fig. [8j where we 
depict the two cases A = ±1. Despite that we can not de- 



FIG. 8: ‘Winding’ of the vector d(k ) for the linearized model 
© for (a) A > 0 and (b) A < 0. 


fine a winding number for neither of the values ±A, we 
can still consider the difference in winding number 5v 
between the two cases, which gives \8v\ = 1. Therefore, 
we expect a zero energy bound state at a boundary be¬ 
tween two regions with A = ±1 respectively, even in the 
linearized model. We stress, that although that this ar¬ 
gument in not rigorous, it is nevertheless true, and in the 
main text we showed that the analytic form of the zero 
mode of the linearized model of Ref. [2dQ accurately de¬ 
scribes the DZM in the junction of the full Kitaev chain. 

We now turn to the alternative linearized model given 
by 1-L V in Eq. ([ 9 ]). Here the fc-space is again not com¬ 
pact, and there is also a discontinuity at k = 0. The 
first issue is remedied by identifying the points at ±00 
(which amounts to considering the a —> 0 limit of the 
lattice model). To deal with the second, we note that 
for this model, d(k) = ^(0, Asgn(fc), — g + vp\k\), and in 
Fig. [9] we show the corresponding ‘winding’ of the vector 
d(k) = d/\d\, in the case A > 0. Even when identify¬ 
ing the points at k = ±00, the curve is not continuous, 
but with a regularization that snroothens out the singu- 
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larity in the V-shaped band, by replacing the factor sgn 
by a continuous odd function that rapidly changes sign 
around k = 0, the d-vector will will be continuous, and 
the winding number will be well defined. This concludes 
the demonstration of the existence of a linearized contin¬ 
uum model with topological properties identical to that 
of the Kitaev chain. 



ua ’w ’=(a ( -V T) • <B3) 

with A(k) = Ho(k) + A. 

Next, we note that Det (H s (k)) = Det(?7 H s (k)W) = 
Det(A(fc))Det(A(— k) T ) so that if H s (k ) is gapped for 
all fc, i.e. Det (H s (k)) ^ 0, the determinant of A(k) 
can not vanish either. This allows us to define z(k) = 
exp (i0(k)) = Det(A(fc))/|Det(A(fc))| for gapped Hamil¬ 
tonians H s {k). One may then show that 


FIG. 9: The ‘winding’ of the d(fc)-vector corresponding to the 
linearized model T-L v before the regularization which removes 
the discontinuity at fe = 0. The arrows indicate the direction 
in which k increases. We have used A > 0. As discussed in 
the text, the gap between the points 0 + and 0~ is closed if 
the dispersion relation is smoothened at k = 0. 


Appendix B: Topological aspects of the s-wave 
paired models 

In this appendix, we discuss the topological proper¬ 
ties of the full and linearized s-wave models. Due to the 
extra spin degree of freedom in these models, the wind¬ 
ing arguments used for the p-wave superconductors are 
not directly applicable and another method of topologi¬ 
cal classification must be used. We will use the method 
outlined in Ref. l22l 

We begin with the full s-wave model, given by Q. 
Assuming a real and constant order parameter, the cor¬ 
responding fc-space Hamiltonian can (in suitable units) 
be written as 


H s (k) = (fc 2 - I-l)t z sq - A TySy (Bl) 

where the Pauli matrices r t and Sj act in particle-hole 
space and spin-space respectively. This Hamiltonian be¬ 
longs to symmetry class BDI, meaning PTRS T 2 = +1 
and PHS C 2 = +1. These operators are in our chosen ba¬ 
sis given by T = tqSqK and C = t x sqK, with K denoting 
the complex conjugation operator. 

To investigate the topological properties of this Hamil¬ 
tonian, we write it in the form 


H s (k) = 


H 0 (k) 


A 


A T -H 0 (k). 


(B2) 


where the matrix structure is in particle-hole space, 
H 0 (k ) = (k 2 — /2)s 0 and A = A is y . We note that the 
latter term is real and has the property A T = — A. By a 
unitary transformation with U = exp(—*(7 t/4)tj,So), the 
matrix in equation (B2) can be rotated into 


z(k) = sgn(A 2 + ( k 2 - nr)- 


(B4) 


For the s-wave Hamiltonian (Bl), which is real and 


gapped for all k, p and finite A, the determinant is real 
and non-vanishing. Then z{k) is well defined and is equal 
to +1, independently of any compactification of fc-space 
(which is needed for any well defined topological invari¬ 
ant), rendering the model topologically trivial. 

We next turn our attention to the linear s-wave model. 
To derive it, we apply the linearization scheme described 
in section |II C to @ and again assume a real order pa¬ 
rameter which gives us the following Hamiltonian: 


#Lin = / dx^2({-iv F ifil + d x ip at+ + iv F tpl_d x (p <7 -)+ 
J CT 

AWOPt.+n,- + + <4,+<4 


We write this as 


with 


njS _ 

^-Lin — 


H lin = J dx^U S Ua {x)^ 


(—ivpd x 0 0 —A ^ 

0 —ivpd x A 0 

0 A ivpd x 0 

y —A 0 0 wf 9 x J 


(B5) 


(B6) 


(B7) 


and the basis 


^ = (V’t.+i ¥4,+> <4-) 1 


(B8) 


The matrix in equation (B7) looks very much like 


two separate blocks of the linear p -wave superconductor 
which seems a bit troublesome since we know that the 
linear p -wave model host zero modes. That would imply 
that the linear s-wave model also would host zero modes, 
which would contradict our findings in this paper. 

One may suspect that the appearance of two p -wave 
models is incidental, and that by adding corrections to 
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the linearization this illusion is shattered. This suspicion 
is indeed justified since, as we now show, the linear s- 
wave superconductor in fact is topologically trivial. 

The corresponding fc-space Hamiltonian is given by 


H hin( k ) = v F kr z s o 


Ar y s y 


(B9) 


where the Pauli matrices r,; and s,; now act in right-left 
space and spin-space respectively. This Hamiltonian also 
belongs to class BDI. In our basis the particle-hole (now 
right-left) symmetry operator is given by C = t z sqK 
and the pseudo time reversal symmetry operator is 7” = 
t x s x K. We rotate the Hamiltonian with the unitary ma¬ 
trix U = exp(— i{ir/£)T x s x ), giving us a structure like 
(B3| but now with A(k) = iv F ks x + iAs y . One may 


then, as above, define z(k) which in this case turns out 
to be z(k) = sgn (y 2 F k 2 + A 2 ) = +1 for all v F , k and 
finite A. As was the case in Appendix [A] k-space is not 
compact. Regardless of this issue, z(k) can never wind. 

Thus we can conclude that both the full and linear 
s-wave superconductors are trivial, and hence that the 
zero modes these models exhibit are not topologically 
protected. 


Appendix C: Origin of the term 


Starting from the original Lagrangian (19), we present 
an argument for how the kinetic term 
can appear in an effective Lagrangian. Although, as al¬ 
ready emphasized in the main text, several of the steps 
in the below derivations are based on unproven assump¬ 
tions, the emergence of the kinetic term is far from obvi¬ 
ous, and this indicates that a more rigorous proof along 
these lines might be possible. 

The starting point is the partition function, 


z[9,g] = / 


i f d 2 x 


(Cl) 


The strategy is to change fermionic variables in such a 
way that the high energy part of the spectrum can still be 
bosonized and integrated out, as in the previous section, 
while the the low lying fermion spectrum will be captured 
by a Lagrangian like (33). To this end, we shall use the 
following identity, 

J V[a y ]V [^, £]e * 1 ^ ’)-««] 

= J T>[a M ]e i ^' l2x “ ,1 - J '' + 5 Trln ( //+a,lp ' i ) 




(C2) 


where p x = —id x , j^ = ipj^ip and U^ is an Hamiltonian 
that we shall assume to be quadratic in the fields and 
H is the corresponding operator acting on the Nambu 
spinors. 


To derive the last line in (C2) we first calculate the 


the trace (which is over both space and Nambu indices). 
The resulting integrals are not convergent in the ultravio¬ 
let since there is no time derivative, so we must introduce 
a cutoff energy scale A. The resulting effective theory is 
only to be applied below this scale. Note that there is 
no gauge invariance related to the auxiliary field a since 
it does not couple to a conserved current. Taking for U^ 
the expression (35) a straight forward calculation gives 
Tr In (H + a^p y ) = coa§ + cia 2 +... where we omitted all 
higher derivative terms. The explicit expressions for the 
coefficients in terms of S, M 2 and A are not particularly 
illuminating. Substituting this in the second line of (C2) 
and integrating over a, we retain the third line with 


F\jn] = Co jo + Cl ii + • ■ • ■ 


(C3) 


Before inserting the identity (C2) in the path integral 


(Cl), we perform the chiral rotation, 


ip ' 


a jl?(x)75 


Ip 


under which the Lagrangian (19) becomes 


C = 1 p(i$~ 7T/ fc - g(x)) 0 • 


(C4) 


(C5) 


Putting this together, we get the following representation 
for the partition function, 

Z[0,g] = , (C6) 

S = J d 2 x [Pp (i$ — F j k — </l — to) ip (C7) 

+ Copipjo'tp) 2 + Ci(?bT0) 2 - - U ? ] , 

where we put g(x) = m to connect to the previous dis¬ 
cussion about the kink solutions. Next we make a shift 
—> a p — nj k , to rewrite the action as 

S = J d 2 x[C f ~ + irj£ - Ht] , (C8) 

where 

Cf = tp (i$ - (/l - m) ip + c 0 (^ 7 oV’) 2 + ci(ipjiip) 2 (C9) 

is very similar to the massive Thirring model. The '0-field 
can now be integrated to give an effective Lagrangian, for 
the field. Using the gauge invariance of (C9) we get 


1 


^e//( a ) — ~ 2 F 


m 


(CIO) 


lowest order by expanding the logarithm and evaluating 


where F llv is the field strength for the potential a y and 
m a dimensional constant that depends both on m and, 
via the coefficients Cq and ci, on S, M 2 and A. Finally, 
we can integrate the vector field a y to get the desired 
effective action for the £-field, 

Cieb = \^d„8 (6 m + tfid^) -Us- b^ + ... (Cll) 
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where we also used the constraint (27) to express j? in 
terms of 0, and where the dots indicate both neglected 
higher derivative terms in the quadratic action, and in¬ 
teraction terms resulting from integrating the field. 
All the steps glossed over above can be performed, at 
least to low order in perturbation theory. The main ques¬ 


tion is however not technical, but rather what principle 
should be used to determine TL j. A possible approach is 
to choose the parameters in so to minimize the size 
of the leading corrections due to higher derivative terms 
and induced interactions. 
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